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E. J. Stone, M.A., F.R.S., President, in the Chair. 

John Wilson Appleton, 12 Amberley Street, Toxteth Park, 
Liverpool; and 

The Hon. Cecil Buncombe, Nawton Grange, Nawton, North 
Yorkshire, 

were balloted for, and duly elected Fellows of the Society. 


On Neivtoris Solution of Kepler's Problem. By Professor 
J. C. Adams, M.A., F.E.S. 

Of all the methods which have been proposed for the solution 
of this problem, that which leads most rapidly to a result having 
any required degree of precision may be briefly explained as 

follows . 

The equation to be solved by successive approximations is 

x - e sin x = 

where z is the known mean anomaly, e the eccentricity, and x 
the eccentric anomaly to be determined. 

Suppose x Q to be an approximate value of ?, found whether by 
estimation, by graphical construction, or by a previous rough 
calculation, and let 

x Q -e sin ~ 0 , 

Then if 
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and 

x f =x 0 + $x 0 . 

x! will be a much more approximate value of 
Similarly, if we put 


A — e sin x J = z\ 

and if 

/ Z Z 

0 X = 

i — e eos x’ 

and 

x" = a?' + 5 x f , 


x lf will be a much more approximate value of x than x r * and so 
on, to any required degree of approximation. 

If the error of the assumed value x 0 be supposed to be of the 
order i, w 7 hen e is taken as a small quantity of the larst order, 
then the error of the value x ’ will be of the order ?.i+ i =i f 
suppose, similarly the error of the value x" will be of the order 
2 i' + i = 4 ^ + 3 , and so on, so that the order of the error is more 
than doubled at each successive approximation. 

The above explains the immense advantage of this process 
over the use of series proceeding according to powers of e, when 
great precision is required in the result; since, in this latter 
method, the addition of a new term only increases the order of 
the error by unity. 

The degree of rapidity of the approximation may be still 
further increased by the following slight modification of the 
above process. 

Starting, as before, with the value x 0f and calling £— z 0 =8z 0 , 
we should obtain a much more accurate value than before of the 
correction 8x 0 to be applied to x OJ by putting 


° l — e cos (x Q + |r $ t37 0 ) i — e cos (x Q + J 5 ^ D ) * 

Now, e being supposed to be small, 8z a is an approximate value 
of Sx of and may be written for it in the small term, in the de¬ 
nominator. 

Hence, if we put 

§ t r _ 5 

l ° I -<? COS (tf 0 -r^- §^ 0 ) ? 


[V ~ x Q + 5 x 0 , 

x f will be a nearer approximation to the true value of x than was 
obtained before by the corresponding operation. 

Similarly, if * 

x f — e sin x f = s ! , 


and 


z — s' 
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and if 


then 


$a/= 


I —e 


cos (x' + £ 5 z')* 


x" ss x f + 5 X f 
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will be the next approximate value of x, and the process may be 
continued as far as we please. 

If the error of x Q be of the order i, that of x l will now be of 
the order that of x" will be of the order 2 (2i + 2)4-2 

= 4 i + 6 , and so on, so that the degree of rapidity of the ap¬ 
proximation is still greater than before. 

If we chose to take the mean anomaly itself as the first 
approximate value of the eccentric anomaly—that is, if we put 


we should have 


z Q = z—e sin z, 

and the value of $x Q given by the first method would be 


B x 


o ~ 


e sin z 
i — e cos z 


while that given by the second and more accurate method would 
be 

5 a? 0 - — -——v—, 

I — e cos (z + £ e sin z) 

and the error of x' = x 0 + 8 x 0 would be of the 3rd order in the 
former case, and of the 4th order in the latter. 

In practice, however, a much nearer first approximate value 
of x may be always found by inspection, and of course the 
smaller the error of this value is, the more rapid will be the rate 
of the subsequent approximations. 

The methods above explained have been long known. The 
first method is given at p. 41 of Thomas Simpson’s “ Essays on 
several Subjects in Speculative and Mixed Mathematics,” pub¬ 
lished in 1740; and Gauss’ method given at pp. 10-12 of the 
44 Theoria Motus,” published in 1809, is essentially the same. 

The second method, or rather the modification of the first, is 
given by Cagnoli in his 44 Trigonometrie,” at pp. 377, 378 of the 
first edition, published in 1786, and at pp. 418-420 of the second 
edition, published in 1808. 

How, my object in the present note is to point out that the 
first method explained above is exactly equivalent to that given by 
Hewtonin the 44 Principia,” at pp. 101, 102 of the second edition, 
and at pp. 109, no of the third edition, when Hewton’s expres¬ 
sions are put into the modern analytical form. 

Hone of the subsequent authors, however, mentions this 
method as being Hewton’s, the unusual form in which Hewton’s 
solution is given having, no doubt, caused them to overlook it. 

F 2 


© Royal Astronomical Society • Provided by the NASA Astrophysics Data System 


Downloaded from http://mnras.oxfordjoumals.org/ at Florida Atlantic University on July 12, 2015 







1882MNRAS. .43. . .43Ai. 


46 Prof. Adams , On Newton’s xliii. 2, 

I11 the first edition of tlie u Principia ” a modification of the 
method is given which was, I have no doubt, intended by 
Newton to be equivalent to the second method given above ; but 
by some inadvertence, instead of the denominator of Saf being 

I — G COS (x r + Jr sy), 

when expressed in the above notation, he takes it to be what is 
equivalent to 

.1 — e cos (x 1 +2 C s ' in x'), 

which is only true for the first approximation when w Q is taken 

In the second and third editions this error is corrected, but 
Newton contents himself with the more simple expression given 
by the first method. 

We need not be surprised that Newton should have employed 
this method of solving the transcendental equation 

x — e sin x — z, 

since the method is identical in principle with his well-known 
method of approximation to the roots of algebraic equations. 

For convenience of calculation, the approximate values 
0? o , x\ x ,f , &c., should be so chosen that their sines may be taken 
directly from the tables without interpolation ; and, since each 
approximation is independent of the preceding ones, this may 
always be done if x' be taken equal, not to 03 o + Sx Q itself, but to 
the angle nearest to x Q + 8 x Q which is contained in the tables, 
and if similarly x" be taken equal to the tabular angle which is 
nearest to x' + 8 x', and so on. In the first approximation it will 
be amply sufficient to use 5-figure logarithms, but in the sub¬ 
sequent ones tables with a larger number of decimal places 
should be employed. 

A first approximate value of the eccentric anomaly corre¬ 
sponding to any given mean anomaly may be found by a very 
simple graphical construction, provided we have traced, once for 
all, a curve in which the ordinates are proportional to the sines 
of the angles represented on any given scale by the abscissae. 

This curve is commonly called “ the curve of sines.” It 
will be sufficient to trace the portion of the curve for which the 
ordinates are positive. 
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Let A 0 B be the line of abscissae, and let A 0 be taken equal 
to 0 B, and let each of them be divided into 90 equal parts 
representing degrees of angle. Let A N be any abscissa re¬ 
presenting the angle x, and let the corresponding ordinate hi P 
= c sin x ; then the greatest ordinate will be OC = c, correspond¬ 
ing to the abscissa AO. 

Suppose the curve line A P C B to be divided into 180 parts 
which correspond to equal divisions on the line of abscissae 
AN OB. 

Then if E be taken in A 0 so that E 0 = e X 57*296 divisions, 
or if AE =90 —6x57-296 divisions, and if CEbe joined and 
PM be drawn parallel to it through P meeting the line of 
abscissae in M, then A M will represent the mean anomaly corre¬ 
sponding to the eccentric anomaly represented by A N. 

For, since the triangles PMN, CEO are similar, 


MN_PN 
EO CO 


= sin 


:c, 


and therefore MN = E 0 sin x =57*296 (e sin as). 

Hence MN represents the number of degrees in as— 2, and 
therefore A M represents the mean anomaly z. 

Conversely, if A M represents any given mean anomaly, then 
if M P be drawn parallel to E C, it -will cut the curve in the 
point P corresponding to the eccentric anomaly. 

By the employment of a parallel ruler we may find the 
eccentric anomaly corresponding to any given mean anomaly, 
or conversely, without actually drawing a line. For if we lay 
an edge of the ruler across the points E C and then make a 
parallel edge to pass through the point M it will cut the curve 
in the point P required. 

Thus we may always find a first approximate value of the 
eccentric anomaly, without making repeated trials, whether the 
eccentricity be large or small. 

I described this graphical method of solving Kepler’s 
problem at the Birmingham meeting of the British Association 
in 1849. It is referred to in a paper by Mr. Proctor in vol. 
xxxiii. of the Monthly Notices , p. 390. 

The construction is so simple that it has probably been pro¬ 
posed before, though I have nowhere met with it. 


Note on Professor ZengePs solution of the same problem given 
in Number 9 of the last volume of the “ Monthly Notices 

The only peculiarity in this solution is in the mode of obtain¬ 
ing the first approximate value employed. The subsequent 
approximations are carried on by means of the first method given 
above. Professor Zenger’s process may be represented in a 
slightly different form as follows 
We have 

x~z~e si n x, 
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